AD-A182  188 


OTIC  FILE  COPJ 

NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


DTIU 


THESIS 

APPLICATION  OF  A  MODIFIED  TIME  DELAY 
SPECTROMETRY  TECHNIQUE  IN  MODELING 
OF  UNDERWATER  ACOUSTIC  PROPAGATION 

by 

Louis  L.  Prudhomme 
March  1987 

Thesis  Advisor:  O.B.  Wilson,  Jr. 

Approved  for  public  release;  distribution  is  unlimited 


UNCLASSIFIED 

•>  C*  Ayi  c  CAT:ON  Of  ^«!  Pag? 


REPORT  DOCUMENTATION  PAGE 

•j  REPORT  SEC-RiTv  COkSSif  CAT, ON 

’b  RESTRICTIVE  VARKiNGS 

UNCLASSIFIED 

.’j  SECUR  Tv  CLASS'f'CAT.ON  AoTmORiTv 

}  DISTRIBUTION/ AVAILABHITV  Of  REPORT 

Approved  for  public  release; 

.’b  DEL  <.  ASS  f  CAT  ON  OOANGRAO'NG  SCHEDULE 

distribution  is  unlimited 

ERfQRVNG  ORGAN  lAT'ON  REPORT  N.jV8£R(S) 


S  VONiTOR'NG  ORGANISATION  REPORT  NuV3E»(S/ 


6b  Off  CE  S'V80l  I  ’*  NAVE  Of  VON* T OR  NG  ORGANiSAT.ON 
lit  tpphaoit)  y 


Code  61 


(k  ADDRESS  Cry  Stitt  tnd  SiPCodt) 


Monterey,  California  93943-5000 


Naval  Postgraduate  School 


"b  ADDRESS  (Cry  Sur*  tnd  JiPCodt) 

Monterey,  California  93943-5000 


3a  NAVE  O'  fuNO'NG*  SPONSOR  NG 

ORGAVSAT  ON 


3c  ADDRESS*  Cry  Stitt  tod  tlf>  Codt) 


8b  Off'CE  S»V80l  9  PROCUREMENT  *NSTRuVENT  OEN' f  CANON  V.V3ER 

(if  ipphabft) 


'0  SOURCE  Of  fi.NOiNG  N'j  VBE  #S 


PROGRAV 
ElEVENT  NO 


PRO.ECT 

rA$* 

NO 

MO 

A'  R«  *,  • 

ACCESSON  NO 


’  *.t  ncludt  Stcurty  OitSificttion) 

APPLICATION  OF  A  MODIFIED  TIME  DELAY  SPECTROMETRY  TECHNIQUE  IN  MODELING 
OF  UNDERWATER  ACOUSTIC  PROPAGATION 


'•a  •  :•  RfROR’ 

Master's  Thesi: 


Mr  •  VE  CO.  E  *E  D 
ffiOV  ’0 


OSA’  COOES 


SwB  G«O.P 


SuSiEC’  TE  RVS  Contmut  on  rtotrtr  if  ntftui'y  tnd  idtn |.Ty  ft,  ft  c<  * 

Multipath  Acoustic  Propagat ion ;  Tine  Delay 
Spectrometry  (TDS^ ;  Surface  Reflection 


9  A3S  ARaC*  'Confmu*  on  rtvtlt  if  ntttittry  ind  identify  by  blot *  nvmbtr) 

The  analysis  of  sound  propagating  by  multiple  pat 
short  ranges  has  been  conducted  using  a  Modified  Time 
(TDS )  technique.  In  this  version  of  TDS ,  a  source  dr 
slide  and  an  HP3561A  Dynamic  Signal  Analyzer  are  used 
amplitude  as  a  function  of  frequency  of  sianals  trave 
paths  and  having  different  arrival  times.  Two  sets  o 
acoustic  test  ranges  at  the  Naval  Undersea  Weapons  Fin 
were  analyzed  for  different  environmental  conditions 
relative  amplitudes  of  the  directly  propagating  and  s 
signals.  Comparisons  with  simple  rough  surface  scatt 
reasonable  agreement.  Results  and  control  software  a 
discussed.  Recommendations  for  future  applications  a 


hs  in  ai 
Delay  * 
iven  by 
to  me  a 
1 i ng  by 
f  data 
g ineer  l 
to  dete 
ur  face 
er  i no  t 
re  pres 
re  made 


Spect 
a  1  l 
sure 
dif  f 
f  ron 
nti  St 
rm  i  no 
ref  *  * 
he  or  y 
en  ted 


san  a t 
:  rome  try 
.near  FV. 
t  he 
:  e  r  e  n  t 
the 
at  ;  *  * r. 

■  t  hu¬ 
rt  ed 
siii  w» 


-  i  •  c*.  a » a  mi'i  o*  abs’rac* 

?£  \  ASS  »  ED'.V  V'ED  □  SAVE  AS  •  »’ 


.  s  NAVE  1  RE  SAG  V,  ScE  ND.O.A, 

Prof.  Q.B.  Wilson,  . 


jr  A81TRAC’  SEC^RC'  -ASS  I 

Unclass l f led 


"?4’Mrw-N?s«r“ 


-,cu  d' 


00  EORM  147),  8«  va» 


8JA - -4»  b« -»«a 

A  c"*’  M  I  o"t  M  cbioiftt 


SK  v»  '  *  i  , 

■  hci.As 


Approved  for  public  release,  distribution  is  unlimited 


Application  of  a  Modified 
Time  Delay  Spectrometry  Technique  in  Modeling 
of  Underwater  Acoustic  Propagation 


by 


Louis  L.  Prudhomme 
Lieutenant,  United  States  Navy 
B. S Tulane  University,  1980 


Submitted  in  Partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ENGINEERING  ACOUSTICS 


from  the 


Author: 


NAVAL  POSTGRADUATE  SCHOOL 
March  1987 


Approved  by: 


B.  Wilson,  Jr. , 


Thesis  Advisor 


Gordon  E.  Schacher 

Dean  of  Science  and  Engineering 


2 


The  analysis  of  sound  propagating  by  multiple  paths  in 
an  ocean  at  short  ranges  has  been  conducted  using  a  Modified 
Time  Delay  Spectrometry  (TDS)  technique.  In  this  version  of 


TDS,  a  source  driven  by  a  linear  FM  slide  and  an  HP3561A 
Dynamic  Signal  Analyzer  are  used  to  measure  the  amplitude  as 
a  function  of  frequency  of  signals  traveling  by  different 
paths  and  having  different  arrival  times.  Two  sets  of  data 
from  the  acoustic  test  ranges  at  the  Naval  Undersea  Weapons 
Engineering  Station  were  analyzed  for  different  environmen¬ 
tal  conditions  to  determine  the  relative  amplitudes  of  the 
directly  propagating  and  surface  reflected  signals.  Compar¬ 
isons  with  simple  rough  surface  scattering  theory  showed 
reasonable  agreement.  Results  and  control  software  are 
presented  and  discussed.  Recommendations  for  future 
applications  are  made. 
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The  reader  is  cautioned  that  computer  programs  developed 
in  this  research  may  not  have  been  exercised  for  all  cases 
of  interest.  While  every  effort  has  been  made,  within  the 
time  available,  to  ensure  that  the  programs  are  free  of 
computational  and  logic  errors,  they  cannot  be  considered 
validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  Naval  Undersea  Warfare  Engineering  Station  (NUWES) 
in  Keyport,  Washington,  performs  measurements  on  torpedo 
radiated  noise  as  part  of  the  test  and  proofing  process  of 
torpedoes.  From  these  measurements  torpedo  radiated  noise 
baselines  are  established,  torpedoes  are  accepted  or 
rejected  on  the  basis  of  meeting  established  radiated  noise 
criteria,  and  design  improvements  in  noise  reduction  can  be 
evaluated.  For  these  reasons,  it  is  important  that  sound 
transmission  loss  between  the  source,  the  torpedo,  and  the 
receiver,  a  system  called  the  Noise  Recording  System  (NRS) , 
be  accurately  modeled.  The  current  model,  based  on 
empirical  data,  appears  to  be  satisfactory  under  normal 
conditions.  However,  under  unusually  rough  surface 
conditions  scattering  effects  may  disturb  the  measurement. 

A  limited  number  of  experiments  and  theoretical  studies 
have  been  conducted  by  Naval  Postgraduate  School  (NPS) 
Students  and  Faculty  to  determine  the  significance  of  rough 
surface  reflection  in  the  multipath  acoustic  transmission 
loss  problem.  Brekke  [Ref.  1]  developed  a  computer- 
controlled  FFT-based  dynamic  signal  analyzer  variation  of  a 
technique  called  Time  Delay  Spectrometry  (TDS)  to  measure 
the  separate  contributions  of  multipath  propagation. 


TDS,  originally  developed  by  Richard  C.  Heyser  in  1967 
[Ref.  2],  utilizes  a  transmitted  Linear  Frequency  Modulated 
(LFM)  pulse  or  "swept  tone"  with  constant  amplitude.  A 
frequency  tracking  spectrum  analyzer  uses  the  difference  in 
arrival  times  to  discriminate  between  the  acoustic  signal 
traveling  by  a  direct  path  and  those  signals  traveling  by 
reflected  paths. 

In  a  thesis  at  NPS,  LT  Ward  [Ref.  3]  outlined  the  multi- 
path  problem  including  environmental  considerations  for  the 
NUWES  test  ranges  at  Dabob  Bay  and  Nanoose.  Further,  he 
generated  a  computer  model  for  predicting  transmission  loss 
under  varying  sea  states.  This  model  has  not  yet  been  veri¬ 
fied  experimentally.  However,  TDS  may  now  make  this 
possible. 

B.  OBJECTIVES 

The  objective  of  this  thesis  is  to  apply  the  modified 
TDS  multipath  measurement  system  to  measure  the  effect  of 
surface  roughness  on  propagation  loss  as  a  function  of 
frequency  and  compare  the  results  with  theory.  The  problem, 
theory,  and  special  considerations  will  be  presented  first, 
followed  by  a  description  of  the  TDS  systems.  Analysis  of 
data  taken  at  the  test  range  at  Nanoose  will  then  be 
presented  including  a  brief  discussion  of  the  variability  of 
the  results.  Results,  conclusions  and  suggestions  for 
future  research  will  also  be  discussed. 


II .  PROBLEM  DESCRIPTION  AND  CONSIDERATIONS 

A.  THE  MULTIPATH  PROBLEM 

For  the  measurement  of  radiated  noise,  during  part  of 
the  test  run,  the  torpedo  passes  by  a  vertical  array  of 
three  calibrated  omnidirectional  hydrophones.  Maximum 
horizontal  range  between  the  torpedo  and  array  is  less  than 
1000  yards.  Torpedo  and  hydrophone  arrays  are  at 
approximately  mid  water  depths.  Water  depth  is  typically 
about  600  feet  in  Dabob  Bay  and  about  1200  feet  at  the 
Nanoose  range.  Transmission  loss  due  to  absorption  can  be 
neglected  since  geometrical  spreading  losses  are  between  40 
and  60  dB  greater  than  absorption  at  these  ranges  at  sound 
frequencies  of  interest  between  50  Hz  and  30  kHz. 

Since  the  tests  are  conducted  in  relatively  shallow 
water,  acoustic  reflection  from  the  boundaries  can  make 
significant  contributions  to  the  total  acoustic  intensity  at 
the  receiver.  Therefore,  in  order  to  model  the  multipath 
problem  adequately,  geometry,  boundary  conditions,  and 
environmental  effects  must  be  considered. 

1 .  Assumptions 

It  is  assumed  that  the  acoustic  signal  is  of  small 
amplitude  and  that  sound  propagation  can  be  described  by  the 
linear  inhomogeneous  wave  equation.  Thus  the  water  channel 
acts  as  a  linear  filter.  It  is  also  assumed  that  the  water 
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column  containing  the  source  and  the  receiver  is 
homogeneous.  The  bottom  is  considered  to  be  flat  and 
horizontal . 

2 .  Multipath  Geometry 

Figure  2.1  depicts  simple  multipath  geometry. 
According  to  the  approach  developed  in  Albers  [Ref.  4:pp. 
49-51]  and  presented  by  Ward  [Ref.  3:pp.  87-88]  and  Brekke 
[Ref.  l:pp.  25-27],  the  acoustic  signal  at  the  receiver  is 
the  sum  of  source  signals  traveling  the  direct  path  and  the 
contributions  of  "image1'  sources  traveling  various  reflected 


paths . 


Following  Brekke's  notation  [Ref.  l:p.  27],  the 


source  depth  is  ZS  and  the  receiver  depth  is  ZR.  The 
horizontal  separation  between  source  and  receiver  is  R  and 
the  channel  depth  is  H.  In  terms  of  these  quantities,  the 
direct  path  distance  (XR) ,  the  surface  reflected  path 
distance  (XSR)  and  bottom  reflected  path  distance  (XBR)  are 
given  by: 


XR  =  [R2  +  (ZR  -  ZS)2]1/2  , 


XSR  =[R2  +  (ZR  +  ZS)2]1/2  , 


(2-1) 


(2-2) 


XBR  =  [R2  +  (2H  -  ZR  -  ZS)2]1/2 


(2-3) 
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Figure  2.1  Simple  Multipath  Geometry 
(Extracted  from  Ref.  1] 
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p.  112].  The  acoustic  pressure  at  the  receiver  due  to  the 
signal  arriving  via  the  direct  path  is  given  by: 

Pd  *  (A/XR)  exp  ] ;  t-(k(XR))l  (2-8) 

The  contribution  of  the  surface  "image"  is  given  by: 

Ps  -  -R,  (A/XSR)  exp  3:  t-(kiXSR)) )  (2-9) 

where  R«  is  the  coefficient  of  reflection  for  the  surface 
Ref.  3 : p .  90  .  Similarly,  the  contribution  due  to  bottom 
reflection  is  given  by: 

Pjj  *  Rjj  (  A  XSR )  exp 3  [  t  -  (  ki  XBR  *  )  ;  1  2  -  1  > , 

where  is  the  coefficient  of  reflection  for  the  bottom. 

It  is  possible  that  have  an  additional  phase  term 

determined  by  the  nature  of  the  water-bottom  composition  and 
the  angle  of  incidence  ;Ref.  3:p.  90 ; .  Assuming  that  the 
direct  and  reflected  waves  are  coherent  one  finds  the  tot  a. 
acoustic  pressure  at  the  receiver  is  the  sum  of  several 
cont  r lbut 10ns : 


ptot  -  pd  ♦  ps  ♦ 

Equation  (2-11)  can  be  rewritten  as: 


ptot  * 


N  Pd 


(2-12 


whera 


I 


I 

3 


X 

sy 


'.V 

,Vi 


Rs  oxp(].ts)  +  ~  Rbexp(].tb) 


(2-13) 


Ward  [Ref.  3:p.  91]  defines  this  factor  as  the  "Multipath 
effect . " 

If  the  phases  of  the  signals  are  incoherent,  the 
intensity  of  the  combined  signals  is  the  sum  of  the 
intensities  of  the  separate  signals,  and  is  given  by: 


tot  *  *d  ♦  ls  >xb 


Equation  (2-14)  can  then  be  rewritten  as: 


(2-14 


1 1 ot  *  MI  *d  - 


(  2  -  1  b  > 


where 


MI  *  t !♦!»/ Id  *  !b/Id] 


(2-16) 


Equation  (2-16)  is  defined  as  the  Intensity  Multipath 
coefficient.  The  surface  reflection  intensity,  Is,  is 
determined  by  the  condition  of  the  surface. 

The  reflection  coefficients,  R^  and  Rb,  are 
determined  by  the  physical  properties  of  the  boundaries 
is  small  at  both  ranges  Refs.  1,  3  and  should  be 
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considered  significant  only  when  either  source  or  receiver 
is  located  near  the  bottom.  For  purposes  of  this  research 
Rb  ^  assumed  to  be  negligible.  The  surface  is  assumed  to 
provide  no  mechanism  for  absorbing  acoustic  energy  and, 
since  the  specific  acoustic  impedance  of  air  is  small 
compared  to  that  of  seawater  [Ref  5:p.  127],  the  energy  in 
an  acoustic  signal  incident  at  the  surface  is  assumed  to  be 
perfectly  reflected  or  scattered  back  into  the  water. 

Urick  [Ref.  6:p.  128]  describes  the  sea  surface  as 

"both  a  reflector  and  a  scatterer  of  sound  and  has  a 
profound  effect  on  propagation  .  .  .  where  source  and 

receiver  lie  at  shallow  depth."  Surface  "roughness," 
defined  in  terms  of  signal  wavelength,  determines  the  nature 
of  the  acoustic  signal  after  interaction  with  the  surface. 
Urick  defines  the  Rayleigh  parameter.  kH  sin  ,  as  the 
criterion  for  acoustic  "roughness,"  where  k  is  the  wave 
number  (  /c) ,  H  is  the  rms  wave  height  (peak  to  trough)  and 
is  the  grazing  angle  [Ref.  6:p.  129].  When  kH  sin  >*  1, 
the  reflection  process  is  called  scattering  and  the  sea 
surface  randomly  scatters  all  the  energy  from  an  incident 
acoustic  signal.  For  kH  sin  <<  1,  the  surface  behaves  as 
an  acoustic  mirror  and  the  reflected  wave  is  completely 
coherent  with  the  incident  wave  [Refs.  3,  6[. 

The  vertical  elevation  of  the  sea  surface  is  a  time- 


variant  random  process  and  is  usually  modeled  statistically. 


A  relatively  good  first  approximation  of  the  probability 
distribution  of  the  surface  wave  height  [Ref.  3:p.  66], 
particularly  when  kH  sin  <  1  [Ref.  7:p.  344],  is  the 
stationary  Gaussian  process.  According  to  developments  of 
Beckmann  and  Spizzichino  [Ref.  8:pp.  73-74],  the  mean  time 
average  scattering  coefficient  <q>  for  a  normally 
distributed  surface  is  given  by: 

<q>  =  X(kz)  q0  (2-17) 

where  qQ  is  unity  (smooth  surface) ,  X(kz)  is  the  spatial 
Fourier  transform  of  the  surface  probability  density 
function  or  "characteristic  function,"  and  kz  is  the  wave 
propagation  vector  for  the  acoustic  signal.  The 
characteristic  function  for  a  Gaussian  distribution  is  given 
by : 


X(kz)  -  exp  (-  2  kz2/2  )  (2-18) 

where  is  the  mean-to-peak  rms  wave  height.  The 
propagation  vector,  kz,  for  the  case  where  the  reflected 
angle  :  r  is  equal  to  the  incident  angle  Jj,  is  given  by: 

:  i 


k 


2k  cos 


(2-19) 


The  term  kz  is  equivalent  to  the  Rayleigh  parameter. 
Beckmann  and  Spizzichino  define  it  as  g1/2  and  use  it 


as  a 


measure  of  surface  roughness  [Ref.  8:p.  82].  Using  the 
definition  of  g  and  Equation  (2-15) ,  the  mean  time  average 
scattering  coefficient  of  Equation  (2-13)  becomes: 

<q>  =  exp(-g/2)  .  (2-20 

The  power  reflection  coefficient,  R^  f  is  the  ratio 
of  acoustic  energy  scattered  by  a  rough  surface  to  the 
energy  scattered  by  a  smooth  surface  (unity)  [Ref.  3:p.  68] 
and  is  defined  by  the  expression: 

R  =  <qq*>  =  exp ( -4k2  j2  cos2  .  (2-21) 

The  surface  amplitude  reflection  coefficient,  Rs,  is  then 
given  by  [Ref.  3:p.  63]: 

Rs  =  R1/2  =  exp(-2k2  -.2cos2  ;^)  .  (2-22) 

Based  on  experimental  determinations,  Clay  and  Medwin  [Ref. 

7:p.  344],  show  that  for  values  of  g*/2/2  greater  than 

unity,  Equation  (2-19)  tends  to  under-estimate  the  true 
value  of  Rs.  They  cite  shadowing  effects  and  variation  in 
actual  wave  shapes  as  reasons  for  these  differences. 
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If  it  assumed  that  there  is  negligible  bottom 


reflection,  the  multipath  term  of  equation  (2-12)  becomes: 

XR 

Hs  =  [i  -  Rg  exp(-ju,ts)  ]  .  (2-23) 

The  sound  pressure  level  at  the  receiver  due  to  sound 
travelling  by  the  direct  path  is  given  by: 


SPLD  =  SPLm  -  10  log10  | Ms | 2 


(2-24) 


where 


SPLm  =  20  log10  |Ptotl  •  (2-25) 

Equation  (2-22)  is  the  total  sound  pressure  level  at  the 
receiver. 

4.  Environmental  Considerations 

Brekke  [Ref.  1]  and  Ward  [Ref.  3]  include  extensive 
environmental  data  for  both  NUWES  test  ranges.  The  key 
environmental  variable  affecting  multipath  propagation  is 
extreme  variation  in  the  sound  velocity  profile  (SVP) . 

There  is  the  possibility  of  strong  negative  gradients 
occurring  above  100  ft  in  depth.  The  resulting  refraction 
may  reduce  delay  time  between  the  signal  traveling  the 
direct  path  and  reflected  paths  so  that  they  may  not  be 
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discernible  [Ref.  l:pp.  27-37]  or  may  affect  the  actual 
angle  of  incidence  at  the  surface. 

B.  TDS  CONSIDERATIONS 

An  analytical  discussion  of  the  modified  TDS  technique 
is  presented  by  Brekke  [Ref.  l:pp.  16-24].  The  points  of 
the  verification  to  be  considered  here  are  the  linear  sweep 
or  LFM  pulse  and  the  response  of  the  HP3561A  Dynamic  Signal 
Analyzer. 

1 .  The  LFM  Pulse 

Figure  2.2  represents  the  LFM  pulse  in  the  time 
domain.  The  amplitude  of  the  signal  is  constant  and  for 
convenience  is  assigned  a  value  of  unity.  The  frequency 
range  F  is  the  difference  between  the  start  frequency  f-^  and 
the  stop  frequency  f2.  The  carrier  frequency  fc  is  defined 
as  the  average  of  f2  and  fj_.  The  sweep  has  a  period  of  T 
seconds  before  it  repeats.  In  terms  of  these  parameters, 
the  instantaneous  voltage  of  the  LFM  pulse  can  be  determined 
as  a  function  of  time  in  the  expression  [Ref.  l:p.  17]: 

X(t)  =  cos[  (TT  (F/T)  t2  -  TTFt  +  2  TTfct]  .  (2-26) 

The  phase  of  the  signal  is: 

=  t(F/T)  t2  +  2  7T  f  xt  . 


(t) 


(2-2 


direct  path  differs  m  frequency  from  the  reflected  signals 
by  an  amount  equal  to  the  sweep  rate,  S,  times  the  delayed 
arrival  time  of  the  reflected  signal  given  in  Equations  (2- 
4)  and  (2-5) . 


2 .  Response  of  the  HP3561A  Dynamic  Signal  Analyzer 

In  the  current  experimental  configuration,  the  LFM 


pulse  is  applied  to  a  broad  band-width  projector.  The 


acoustic  wave  propagates  by  one  or  more  paths  through  the 
water  to  a  hydrophone.  The  electrical  signal  from  the 
hydrophone  is  recorded  in  analog  form  on  magnetic  tape.  Th 
recorded  signal  is  subsequently  analyzed  using  the  HP3561 
Dynamic  Signal  Analyzer.  The  measurements  are  based  on  a 
1024  point  FFT  algorithm.  The  displayed  spectrum  contains 
400  bins.  Frequency  resolution  of  the  analyzer  is 
determined  by  the  selected  frequency  span  F  of  the  analyzer 
in  terms  of  the  bin  separation  b: 


b  =  F/400  Hz  .  (2-30 

For  a  frequency  span  of  5  kHz,  the  corresponding  frequency 
resolution  is  12.5  Hz.  The  time  record  length,  the 
observation  time  tQks  during  which  the  analyzer  receives 
data  from  the  signal,  is  given  by: 


fcobs  =  Vb  sec  . 


(2-31 


This  is  the  minimum  amount  of  time  required  by  the  analyzer 
to  obtain  sufficient  samples  to  carry  out  the  FFT  algorithm. 
The  energy  contained  in  each  of  the  400  bins  is  calculated 
and  made  available  for  display  by  the  analyzer  during  the 
observation  time.  In  most  operating  modes,  the  HP3561A 
displays  only  the  most  recent  time  record  data.  In  Time 
Capture  mode,  however,  the  analyzer  acquires  and  stores  data 
contained  in  40  contiguous  time  records.  Time  Capture  will 
be  discussed  in  more  detail  in  the  next  chapter. 

The  consequence  of  these  relations  is  that  the  delayed 
arrival  times  given  by  Equations  (2-4)  and  (2-5)  must  be 
greater  than  tobs  in  order  to  resolve  the  directly 
propagating  and  reflected  signals  [Ref.  l:p.  48].  For  delay 
arrival  times  less  than  tobs,  the  difference  between  the 
instantaneous  frequencies  of  the  direct  and  reflected 
signals  is  less  than  the  resolution  of  the  analyzer  and 
thus,  the  energy  of  both  direct  and  reflected  signals  are 
contained  in  the  same  bins. 

The  possibility  of  energy  leaking  into  an  adjacent 
bin  and  creating  additional  measurement  errors  is  reduced  by 
using  a  "window  function"  feature  on  the  HP3561A.  The 
result  of  using  a  window  function  is  a  minimum  bandwidth 
separation,  f_band,  given  by  [Ref.  l:p.  48]: 


f  band  =  C  x  F 


(2-32) 


where  C,  a  constant  determined  by  the  window  function 
selected,  is  expressed  in  terms  of  the  frequency  span.  For 
a  Hanning  window,  C  is  0.00375.  Thus,  if  the  selected 
frequency  is  5  kHz,  the  resulting  f_band  is  18.75  Hz.  So, 
in  order  to  resolve  the  direct  and  reflected  signals,  it  is 
necessary  to  insure  that  the  delayed  arrival  times  are 
sufficient  to  maintain  the  minimum  bandwidth  separation. 

This  is  accomplished  by  setting  a  lower  bound  on  the  sweep 
rate,  Sm,  given  by  [Ref.  l:p.  49]: 

Sm  >  f_band/t_delay  ,  (2-33) 

where  t_delay  is  given  by  Equations  (2-4)  through  (2-6) . 

For  surface  reflection,  t_delay  =  ts. 
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III.  TPS  MULTIPATH  MEASUREMENT  SYSTEM 


The  TPS  Multipath  Measurement  System  is  comprised  of  two 
major  subsystems  [Ref.  l:p.  39]: 

(1)  Range  and  test  configuration 

(2)  Pata  processing  system. 

The  range  and  test  configuration  described  here  is  specific 
to  the  data  used  in  this  research.  The  test  format  was 
developed  by  John  G.  Burwell  of  the  Acoustics  Pivision, 

NUWES ,  Keyport,  and  has  been  implemented  by  Keyport 
personnel  in  previous  TPS  research.  The  data  processing 
system  presented  by  Brekke  [Ref.  l:pp.  42-61]  has  been 
modified.  A  description  of  the  hardware  configuration  is 
provided  here.  The  software  used  is  a  slightly  modified 
version  of  the  TPS  software  documented  by  Brekke  [Ref.  1]. 
The  key  features  of  the  TPS  programs  will  be  discussed  here 
and  software  modifications  will  be  documented. 

A.  RANGE  ANP  TEST  CONFIGURATION 

The  data  used  in  this  study  were  obtained  at  the  NUWES 
range  at  Nanoose,  Canada,  in  February  1986.  The  acoustic 
projector  for  the  test  was  the  Sonar  Acoustic  Test  System 
(SATS) .  The  receiver  was  the  Noise  Recording  System  (NRS) 
with  three  calibrated  hydrophones  at  depths  of  330  ft. ,  390 
ft.,  and  430  ft.  Horizontal  separation  between  SATS  and  NRS 


was  approximately  300  yds.  Wave  heights  were  estimated  at  1 
to  2  ft  peak-to-peak  at  the  time  of  the  measurements. 

The  test  signal  was  a  linear  sweep  from  50  Hz  to  30  kHz. 
Three  tests  were  conducted  with  sweep  periods  of  1  second,  2 
seconds  and  4  seconds.  Thirty  sweeps  were  transmitted  for 
each  test. 

The  sound  signals  arriving  at  each  of  the  NRS 
hydrophones  were  recorded  on  separate  channels  of  a  magnetic 
tape  recorder. 

B.  THE  DATA  PROCESSING  SYSTEM 

1 •  Hardware 

The  equipment  configuration  used  is  shown  in  Figure 
3.1a.  It  is  essentially  the  same  as  that  used  by  Brekke 
[Ref.  l:p.  42-43],  with  an  HP  thinkjet  printer  replacing  the 
HP2671G  printer.  For  the  routine  tests,  the  graphics  and 
printing  requirements  could  be  satisfied  using  the  HP 
thinkjet  as  shown  in  Figure  3.1b. 

The  measurements  of  the  signal  are  carried  out  by 
the  HP3561A  Dynamic  Signal  Analyzer  on  the  recorded  signal 
under  computer  control.  Previously,  triggering  the  Dynamic 
Signal  Analyzer  was  done  using  a  1  pulse-per-second  signal 
derived  from  the  output  of  the  time  code  generator  [Ref.  1]. 
This  was  not  possible,  since  the  period  of  the  LFM  pulses 
for  two  of  the  tests  exceed  1  second.  To  accommodate  the 
longer  sweep  periods,  the  trigger  signal  of  an  HP3314A 


Arbitrary  Function  Generator  was  set  to  •:  sec  or.  :  ;  ir.  i 
recorded  on  the  edge  track  of  the  data  tape. 

2 .  Software 

Brekke  [Ref.  1]  provides  a  detailed  documentation 
the  TDS  1  program.  Here,  a  general  discussion  of  the 
capabilities  of  the  TDS  software  will  be  presented  with 
emphasis  placed  on  the  features  used  in  this  research.  The 
modifications  used  in  the  current  version  3C  of  the  TDS  l 


The  signal  is  analyzed  by  one  of  three  methods: 

Quick  Analysis,  Normal  Analysis,  and  Auto  Analysis. 

Referring  to  Figure  3.2,  one  sees  that  these  three  methods 
form  the  three  major  subprograms  of  TDS  1 .  These  analysis 
methods  are  discussed  in  the  following  paragraphs, 
a.  Quick_Analysis 

In  this  method  the  HP3561A  is  operated  in  the 
"Time  Capture  Mode".  The  analysis  is  based  on  data  taken 
from  40  time  records  stored  in  the  time  buffer  of  the 
analyzer  memory.  Each  time  record  consists  of  the  data 
processed  during  one  observation  time  t^g.  The  total  time 
during  which  the  analyzer  processes  data  is  40  x  tobs. 

Figure  3.3  illustrates  the  storage  of  the  LFM  pulse  in  the 
time  buffer.  The  data  are  displayed  as  a  3-dimensional 
magnitude  map.  Each  map  consists  of  60  traces.  The  overlap 
(or  increment)  processing  selected  determines  the  time 
increment  represented  by  each  subsequent  trace.  For  the  20% 
overlap  processing  used  in  the  TDS  program,  each  trace 
represents  a  0.2  x  t0bS  increment  of  the  delay  time.  The 
analyzer  uses  20%  new  data  taken  from  the  subsequent  time 
record  for  each  new  trace. 

The  purpose  of  the  overlap  processing  is  to 
optimize  the  useful  information  displayed  by  the  analyzer. 

In  selecting  an  appropriate  overlap,  the  sweep  rate,  the 
analyzer  frequency  span,  and  the  frequency  span  of  the  LFM 
pulse  must  be  considered.  For  example,  the  4  second  LFM 
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i ;  re  ct  mi  ret.  e  ted  siqnals  is  approximately  r*0  which 
ex.eeis  the  <•••  traces  that  can  be  displayed  by  the  analyzer. 

,u:ck  Analysis  provides  the  operator  a  qu :  k 
overview  of  a  "Spectral  Segment"  ot  a  sinqle  sweep  ot  the 
signal  to  be  analyzed.  A  "spectral  segment"  is  the  portion 
'f  the  I FM  pulse  that  can  be  displayed  by  the  analyzer.  If 
the  start  frequency  of  the  analyzer  is  500  Hz  and  the 
frequency  span  selected  is  5  kHz,  then  only  the  "segment"  of 
the  L.FM  pulse  between  these  two  frequencies  can  be 


processed . 


From  the  spectral  1-D  map,  the  arrival  times  of 


the  "Spectral  Segment"  for  both  direct  and  reflected  siqnals 
van  be  determined  by  adqustinq  the  start  time  of  the 
analyzer.  The  time  when  the  reflected  signal  vanishes  from 
the  "Spectral  Segment"  should  also  be  determined.  These 
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values  are  used  as  input  parameters  for  Normal  Analysis  and 
Auto  Analysis.  After  the  determination  of  the  arrival  times 
and  "upper  limits  of  the  delayed  start  time",  or  time  when 
the  reflected  signal  vanishes,  the  program  displays 
individually  all  magnitude  maps  between  these  time  limits 
for  operator  input.  The  operator  selects  peaks  for  direct 
mi  reflected  signal  using  the  cursor  on  the  Spectrum 
Analyzer.  These  values  are  recorded  by  the  program  and  are 
printed  out  at  completion  of  the  analysis. 

It  should  be  noted  here  that  Quick  Analysis 
requires  only  a  single  trigger  since  only  one  sweep  is 
captured.  The  requirement  for  a  sequence  of  synchronous 
triggers  exists  when  Normal  Analysis  or  Auto  Analysis  will 
be  performed. 

b.  Normal  Analysis 

In  this  method  8  sweeps  are  RMS  averaged  to 
reduce  effects  of  signal  fluctuation  and  increase  signal  to 
noise  ratio  before  making  measurements.  The  measurements 
are  made  in  a  manner  similar  to  those  in  Quick  Analysis  in 
which  the  operator  must  select  peak  values  using  the  cursor 
on  the  HP3561A.  The  program  increments  the  delay  step  after 
each  measurement  to  measure  the  frequency  bin  of  the  later 
arrival  time.  These  measurements  are  repeated  for  every 
delay  step  in  the  time  interval  determined  by  the 
Quick  analysis  method. 
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The  delaystep  should  increment  each  measurement 
by  a  frequency  larger  than  the  frequency  resolution  of  the 
analyzer . 

c.  Autoanalysis 

This  method  automates  the  measurements  performed 
by  NormalAnalysis  method.  Due  to  a  high  noise  level  at  the 
lower  end  of  the  spectrum,  the  operator's  observation  was 
necessary  to  discriminate  between  direct  and  reflected 
signals.  Consequently,  the  AutoAnalysis  method  was  not 
used  in  this  study. 

d.  Analysis  Products 

Each  of  the  methods  produces  a  table  of  voltage 
levels  for  direct  and  reflected  signals  as  well  as  a  noise 
level.  The  program  enables  the  operator  to  plot  the 
magnitude  of  direct  and  reflected  signals  as  a  function  of 
frequency.  This  requires  that  all  the  tabulated  frequency 
values  be  monotonically  increasing.  This  is  not  always  the 
case,  since  often  the  reflected  path  peak  is  recorded  at  a 
lower  frequency  than  the  previous  measurement  due  to  the 
uneven  nature  of  boundary  reflection.  Using  a  longer  delay 
time  will  reduce  this  problem,  but  with  the  result  that 
fewer  samples  are  recorded. 


m 

IV.  DATA  ANALYSIS 


The  analysis  of  a  "spectral  segment"  of  LFM  pulse  with 
the  4  second  period  is  presented  here.  Procedures  used  here 
are  the  same  for  the  analysis  of  any  "spectral  segment"  of 
any  LFM  pulse.  Multipath  problem  considerations  are 
discussed  with  respect  to  TDS  parameters.  Measurements  made 
using  the  Normal_Analysis  method  are  used  to  determine  the 
effects  that  surface  reflection  or  scattering  have  on  the 
acoustic  field  at  the  receiver.  Measurements  by  the 
Quick_Analysis  method  on  several  different  sweeps  are 
compared  to  the  Normal_Analysis  results  to  estimate  the 
sweep  to  sweep  variability  in  results. 

A.  ANALYSIS  PARAMETERS 
1 .  Time  Delays 

The  sound  velocity  profile  given  in  Figure  4.1  is 
typical  for  the  Nanoose  range  in  February  [Ref.  9].  Because 
the  effects  of  refraction  due  to  this  profile  are  small  for 
horizontal  ranges  of  a  few  hundred  yards,  the  water  column 
will  be  assumed  to  have  a  constant  sound  speed  of  c  =  4850 
ft/sec.  The  multipath  parameters  for  the  problem  are  given 
in  Table  1. 

The  minimum  time  delay  given  by  Table  1  is  82  ms. 
The  minimum  frequency  span  setting  for  the  HP3561A  is  5  kHz 
corresponding  to  t0bS  =  80ms. 
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Figure  4.1:  Sound  Velocity  Profile  (Nanoose— February ) 
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TABLE  1 


DISTANCES  AND  TIME  DELAYS 


RCVR  (depth)  XR 

upper  (330)  100 

middle  (390)  104 

lower  (430)  210 


(yds) 

XSR 

(Yds) 

yds 

233 

yds 

yds 

251 

yds 

yds 

266 

yds 

Source  depth:  300  ft 

Horizontal  Range:  100  yds  (300  ft) 

c  =  4850  ft/sec. 


ts  (ms) 
8  2  ms 
91  ms 
96.4  ms 


2 .  The  Rayleigh  Reflection  Coefficient 

Assuming  transmission  loss  to  be  due  to  only 
spherical  divergence,  the  sound  pressure  level  at  the 
receiver  due  to  the  direct  signal  is  given  by: 

i 

i 

i 

SPLS  =  SPL [ P ( 1) ]  -  20  log  XR  .  (4-1) 

The  sound  pressure  level  due  to  a  surface  reflected  signal 
is  assumed  to  be  given  by: 


SPLR  =  SPL  [P(l)]  -  20  log  XSR  +  10  log  (Rs)  .  (4-2) 

P(l)  is  the  sound  pressure  at  unit  distance  from  the  source, 
subtracting  (4-2)  from  (4-1)  yields: 


10  log  Rs  =  SPLd-SPLr  +  20  log  (XSR/XR)  . 


(4-3) 


If  all  terms  on  the  right  hand  side  of  (4-3)  are 
known,  then  a  value  for  Rs  can  be  determined.  Using 
equation  (2-21),  one  can  solve  for  the  mean-to-peak  wave 
height.  For  six  sweeps  analyzed  using  the  Quickanalysis 
method,  these  calculations  gave  an  average  value  for  o  of 
0.54  ft  at  frequencies  less  than  1200  Hz.  The  standard 
deviation  at  these  frequencies  was  on  the  order  of  0.05  ft. 
At  frequencies  above  1200  Hz,  the  calculated  values  for  z 
decreased. 

These  results  are  consistent  with  data  given  by  Clay 
and  Medwin  [Ref.  7:p.  344].  Since  the  wave  heights  were 
estimated  between  1  and  2  ft,  0.54  feet  for  o  seems  to  be 
reasonable.  So  for  the  data  used  here,  the  measurements  made 
by  the  Modified  TDS  Measurement  System  appear  to  be  valid. 

3 .  Selecting  Analyzer  Settings 

Since  the  minimum  time  delay  is  82  ms  (Table  1) ,  the 
maximum  observation  time  to^s  which  permits  the  HP3561A  to 
resolve  the  direct  and  reflected  signals  is  80  ms.  This 
corresponds  to  a  5-kHz  spans  on  the  analyzer.  As  pointed 
out  in  Chapter  II,  choosing  the  maximum  observation  time 
gives  the  best  frequency  resolution  on  the  analyzer. 

Figure  4.2  is  the  3-D  map  output  of  Quick_Anal ysis . 
The  start  frequency  selected  here  is  500  Hz.  Ideally,  each 
5-kHz  span  in  the  sweep  would  be  examined  successively. 
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Figure  4.2:  Quick_Analysis  3-D  Spectral  Map 


B.  MEASUREMENTS 


The  received  Acoustic  Level  RL(f)  at  the  face  of  the 
hydrophone  is  given  by  [Ref,  l:p.  74]: 

RL ( f )  =  ML ( f ) -CAL (f)-HS(f)+IL(f)-47 -AG 


where : 


ML(f)  is  measured  level  in  (dB  re  IV) 

CAL(f)  is  measured  level  of  the  calibration  signal 
HS(f)  is  hydrophone  sensitivity 
IL(f)  is  calibration  insertion  loss 
AG  is  the  Auto  Gain  of  the  Amplifier. 

The  levels  ML( f )  for  these  measurements  are  provided 
in  outputs  from  the  TDS1  program  in  Appendix  B.  For 
selected  frequencies  used  in  this  thesis,  HS(f)  and  IL(f) 
are  found  in  Table  2.  These  values  are  from  the  calibration 
records  for  the  equipment  provided  by  NUWES.  The  values  for 
AG  and  CAL(f)  used  here  are  +40dB  and  OdB.  Since  kH  sino  is 
approximately  unity  at  1.2  kHz,  some  assumptions  about 
signal  reflection  at  the  surface  can  be  made.  For  f  <  1.2 
kHz,  the  effect  of  specular  reflection  from  the  surface  is 
determined  by  Equation  (2-22).  The  upper  and  lower  bounds 
of  the  acoustic  pressure  at  the  receiver  are  given  by: 
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-  2  0 . 1 

-18  1. 

4  kHz 

-18  2.4 

—  0  0 

.0  -182.4 

-20. 2 

-18  2.  • 

5  kHz 

-  132  .  3 

-2  0 

.1  -182.0 

-  2  o  .  2 

-13  2..- 

*  Data 

provided 

by 

Nl’WES 

p  wax  -  PD<1  *  virp  F's  ■■■i-h 

and 

pmin  =  pD^f  ~  v^p  ps'  '  (•*“■' 

For  f  >  1.2  kHz,  the  effect  of  scattering  and 

incoherent  addition  from  the  surface  on  the  average  acoustic 
intensity  at  the  receiver  is  given  by  Equation  (2-15)  where: 


Mj-  'l  +  Ig/Ijjj.  (4-5) 

The  ratio  Is/Ir  is  equal  to  the  ratio  of  measured 


levels  of  both  direct  and  reflected  signal. 


Since  the 


where  L'.’L  is  the  measured  level  in  a  given  frequency  band 
and  N  is  the  number  of  measurements.  Table  4  gives  the  mean 
level  <X> ,  and  "normalized"  level  of  the  standard  deviation, 
10  log(S)  -  <X> ,  for  several  frequencies. 

The  Rayleigh  reflection  parameter  varies  from  about 
1  to  5.5  over  the  frequency  span  in  the  experiment.  The 
energy  averaged  Quick_Analysis  measurements  are  close  to  the 
measurements  obtained  by  the  Normal_Analysis  method.  This 
indicates  that  the  Normal_Analysis  measurements  provide  a 
fairly  good  approximation  to  the  average  acoustic  energy  in 
each  band  of  the  measured  signal. 


TABLE  4 

STATISTICAL  COMPARISON  OF  SURFACE  REFLECTED  LEVELS 


Freq  (Hz) 

Quick  Analysis 

<X> 

(dB  re  IV) 
Normal  Analysis 

101og10 (s) 

1000 

-12 . 6 

-12.0 

-13.4 

1500 

-13 . 9 

-14 . 6 

-  3.8 

2000 

-18 . 4 

-18 . 3 

-  0.25 

2500 

-25.7 

-26. 1 

-  2.4 

3000 

-21.7 

-21.9 

-  1.5 

4000 

-24 . 0 

-21 . 5 

-4.0 

5000 

-24 . 4 

-22 . 9 

-  2.7 
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V.  DISCUSSION  AND  CONCLUSIONS 


A.  DISCUSSION 

1 .  Procedure 

Two  factors  which  caused  the  greatest  difficulty  in 
analyzing  data  were  the  way  in  which  the  trigger  was  set  and 
the  choice  of  the  number  of  sweeps  for  each  data  run.  As 
discussed  earlier,  the  Dynamic  Signal  Analyzer  was  triggered 
by  signals  recorded  on  the  edge  track  of  the  data  tape.  The 
Honeywell  5600E  magnetic  tape  recorder  has  a  total  of  14 
channels.  This  includes  13  data  channels  and  one  edge 
channel,  typically  used  for  voice  comments.  The  edge 
channels  can  be  recorded  without  disturbing  the  data 
recorded  on  the  data  channels.  The  edge  tracks,  located  on 
the  "edge"  of  the  tape,  sometimes  suffer  degradation  caused 
by  tape  guides  on  the  recorder.  The  trigger  for  the  HP3314 
Arbitrary  Function  Generator  was  set  to  4  seconds  and 
recorded  on  the  edge  track.  Several  attempts  were  made  to 
synchronize  the  trigger  with  the  start  of  the  LFM  pulse. 

The  data  analyzed  in  this  thesis  required  approximately  1.9 
second  trigger  delay.  For  the  overlap  processing  used  in 
the  TDS  program,  only  80  percent  of  the  LFM  pulse  was 
processed  in  Quick  Analysis.  Since  the  trigger  was  not 
synchronized  to  the  start  of  the  LFM  pulse,  20  percent  of 
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the  s:  ir.  1 1  ketveen  3  kHz  and  12  kHz  could  not  be  captured  i  y 
the  analyzer. 

The  trigger  level  remained  high  for  approximately 
1.7  seconds.  This  made  Quick.  Analysis  possible  for  parts  of 
the  sweep  arriving  during  times  of  high  trigger. 

Normal  Analysis  was  not  possible,  however,  because  the 
analyzer  continued  to  trigger  and  sample  di f ferent  portions 
of  the  LFM  pulse. 

For  further  TDS  measurements  a  sharp,  short-duration 
trigger  is  necessary  and  should  be  recorded  on  tape  with  the 
data.  The  data  Brekke  used  had  a  one  second  sweep  period 
which  was  synchronous  with  the  range  time  signals.  A  1-Hz 
trigger  was  derived  from  the  IRIG-B,  time  code  generator 
translator.  If  a  one  second  sweep  period  is  maintained,  the 
frequency  span  of  the  signal  will  have  to  be  reduced  to 
accommodate  a  slower  sweep  rate.  If  a  slower  sweep  is  used, 
such  as  the  4  second  sweep  used  here,  the  advantage  of  a 
larger  frequency  span  is  gained.  In  this  case,  a  trigger 
signal  synchronous  with  the  start  of  the  LFM  pulse  should  be 
recorded . 

The  total  data  time  fur  all  three  tests  was  3.5 
minutes,  with  two  minutes  the  maximum  total  time  for  slowest 
sweep  rates.  NormalAnalysis  requires  8  sweeps  to  make  one 
measurement.  For  thirty  sweeps,  a  maximum  of  3  measurements 
can  be  made  before  the  tape  must  be  rewound,  repositioned 


and  started  again.  For  the  5-kHz  span  used  here 


approximately  42  measurements  are  required  and  consequently 
the  rewind  process  was  repeated  a  minimum  of  14  times.  Wit 
five  minutes  of  data  recorded,  the  rewind  process  would  be 
repeated  only  5  times.  Since  the  rewind  process  totals 
about  3  minutes  the  analysis  time  would  be  reduced 
significantly.  An  additional  benefit  is  that  increased 
information  available  from  the  longer  record  could  be  used 
to  determine  the  variability  of  the  surface  reflection.  In 
future  tests,  the  data  recorded  should  be  of  a  longer  time 
at  one  sweep  rate  with  5  minutes  of  data  accumulation  as  a 
minimum. 

2 .  Software 

The  modifications  made  to  the  TDS1  program  have 
made  it  more  flexible  in  making  analysis  of  data. 
Eliminating  a  requirement  for  the  plotter  shown  in  Figure 
3.1  (b) ,  would  streamline  the  process  considerably. 

Although  HP  Thinkjet  graphic  printouts  are  not  as  detailed 
or  colorful  as  the  plots  of  the  HP7470A  plotter,  for 
perishable  measurements  the  HP  Thinkjet  products  are 
adequate . 

The  QuickAnalysis  printouts  given  in  Appendix  B.l 
read  in  4  ms  increments;  this  corresponds  to  16  ms  of  time 
in  the  NormalAnalysis  printouts,  Appendix  B.2. 

The  program  in  current  form  can  be  used  for  routine 
measurements.  Minor  modifications  need  to  be  implemented 
for  the  changes  just  discussed,  but  personnel  who  would  use 


the  products  and  conduct  the  measurements  should  determine 
the  final  form  of  the  program. 

B.  CONCLUSIONS 

The  Modified  TDS  Measurement  System  presented  by  Brekke 
[Ref.  1]  provides  useful  information  about  an  acoustic 
signal  propagating  by  direct  and  surface  reflected  paths. 
When  the  sea  surface  has  some  degree  of  roughness, 
measurements  by  this  method  permit  the  effect  of  the  surface 
reflected  signal  and  the  direct  signal  to  be  considered 
separately.  Further  modifications  to  the  measurement  system 
are  still  required.  In  its  current  configuration,  the 
system  is  cumbersome.  The  removal  of  the  requirement  for  a 
plotter  could  streamline  the  analysis  for  routine  tests. 

The  HP3561A  Dynamic  Signal  Analyzer  can  perform  an 
energy  average  of  a  recorded  bandwidth.  This  technique  may 
be  more  accurate  than  the  current  method  of  recording  only 
the  peak  value,  particularly  for  QuickAnalysis.  In 
NormalAnalysis  the  averaging  over  8  sweeps,  eliminates  most 
of  the  variability  so  there  should  not  be  a  significant 
difference  between  the  average  of  peaks  and  the  energy 
average.  In  Quick_Analysis ,  however,  the  averaging  over  the 
bandwidth  would  improve  the  accuracy  of  measurements, 
particularly  for  the  reflected  signal  at  higher  frequencies. 
The  results  obtained  here  are  not  conclusive.  More  data 
needs  to  be  analyzed.  The  results  reported  here  were  made 
for  a  limited  frequency  range.  Results  for  different  sweep 
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rates  and  frequency  spans  should  be  investigated  and 
compared. 

In  any  future  tests,  more  extensive  measurements  and 
recording  of  environmental  conditions  need  to  be  made.  If 
reflected  signals  are  distinguishable  and  if  more  than  one 
is  present,  the  operator  may  use  TDS_1  to  compare  any  of  the 
reflected  signals  with  the  directly  propagating  signal. 
Future  modification  to  the  TDS  program  should  allow  for  the 
analysis  of  more  than  one  reflected  signal. 


APPENDIX  A 


RANGE  AND  TEST  DATA 

This  list  provides  a  guideline  of  information  to  be 
recorded  at  the  range  just  prior  to  conducting  tests  using 
the  Modified  TDS  Measurement  System. 

1.  ENVIROMENTAL 

(a)  Wind-Speed  and  Direction 

(b)  Sound  Velocity  Profile 

(c)  Wave  Heights 

2.  CONFIGURATION  OF  TEST 

(a)  Horizontal  Ranges  Between  Source  and  Receivers 

(b)  Depths  of  Source  and  Receiver 

(c)  Range  to  Nearest  Land 

(d)  Water  Depth  at  Test  Sight 

3.  TEST  PARAMETERS 

(a)  Calibration  Signal 

(b)  Auto  Gain  Settings 

(c)  LFM  Pulse 

-  Frequency  Span 

-  Repetition  Period 

-  Number  of  Sweeps 

4.  CALIBRATION  RECORDS  FOR  RECEIVING  EQUIPMENT 

(a)  Hydrophone  Sensitivities 

(b)  Insertion  Loss  Data 
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EPRODUCED  AT  GOVERNMENT  EXTENSE 


NPS 


MIDDLE  HYDROPHONE 
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ANALYIS  METHOD :  QUICK  ANALYSIS 

DATE  OF  ANALYSIS:  11/30/86 


RECORDED  MAGNITODE  LEVELS  ( dB  re  IV) 


DELAY(  K S  ) 

DIRECT  KAGN 

DIRECT  PRQ(  HZ  ) 

REFL  KAGN 

REFL  PRQ{ HZ ) 

NOISE  MAGW 

1930 

'6.60 

563 

19.75 

0 

-47 .81 

1934 

-7.00 

575 

19.75 

0 

-46 . 91 

1938 

'4.90 

750 

19.75 

0 

-46 . 20 

1942 

-4.90 

775 

19.75 

0 

-45.30 

1 946 

-6.52 

925 

19.75 

0 

-43 . 57 

1950 

-5.85 

1075 

19.75 

0 

-41.11 

1954 

-4.62 

1213 

-10. 46 

537 

-38 . 72 

1958 

-2.35 

1275 

-  11.74 

613 

-36.62 

1962 

-2.17 

1350 

-12.51 

775 

-36.36 

1966 

-3.96 

1400 

-8 . 80 

838 

-37.98 

1970 

-7.69 

1575 

-10.68 

912 

-39.35 

1974 

-7.91 

1713 

-9. 84 

1075 

-38.47 

1978 

-  7.56 

1775 

-7.72 

1130 

-37.48 

1982 

•7.93 

1250 

-7.03 

1250 

-37.80 

1986 

-  11.70 

1950 

-8.40 

1313 

-39.74 

1990 

10.60 

206  3 

-12.51 

1363 

-41.76 

1994 

-  12.05 

2125 

-19.36 

1525 

-42.98 

1998 

14.31 

2225 

-  15.91 

1713 

-42.02 

2002 

18.56 

2375 

-  14.30 

1800 

-42.52 

2006 

17.16 

2563 

-  11.77 

1913 

-41.81 

2010 

14.03 

2625 

-10. 19 

1975 

-40. 39 

2014 

-  11.78 

2712 

-12.24 

2075 

-40.37 

20  18 

-9.89 

2800 

-10.73 

2162 

-41.09 

2022 

-9  97 

2888 

-18.77 

2263 

-42 . 49 

2026 

11.22 

2975 

-26  56 

2388 

-42.87 

2030 

12.62 

3075 

-21.96 

2513 

-42 . 76 

2034 

14.96 

3250 

-20.09 

2600 

-42 . 19 

2038 

-12.72 

3325 

-20.20 

2700 

-41.66 

2042 

-12.21 

3430 

-17.33 

2787 

-40.95 

2046 

-12.87 

3480 

-21.00 

2838 

-40.82 

2050 

-14.32 

3625 

-20.61 

3075 

-40.81 

2054 

-13.76 

3738 

-16.74 

3100 

-40.84 

2058 

-13.21 

3825 

-18.80 

3238 

-40.11 

2062  - 

14.  37 

3930 

-13.67 

3337 

-39.70 

2066 

13.54 

4025 

-  16.05 

3400 

-39.76 

2070 

12.61 

4150 

-22.63 

3475 

-39.78 

207  4 

12.24 

4213 

-17.31 

3663 

-39.95 

2078 

11.85 

4350 

-  16.63 

3738 

-38.92 

2082 

-  12.54 

4400 

-15.81 

3825 

-30.12 

2086 

-12.14 

4550 

-  16.00 

3938 

-37.79 

2090 

-  12.81 

4625 

-17.19 

4088 

-37.49 

2094 

-  12.08 

4750 

-  14.96 

4188 

-37.41 

2090 

13.14 

4850 

-13.11 

4250 

-38 . 29 

2102 

13.76 

4950 

-13.13 

4350 

-38 . 96 

2106 

-  13.53 

5038 

-  13.67 

4400 

-38.23 

2110 

-  12.81 

5200 

-15.56 

4550 

-38.56 

2114 

-12.17 

5225 

-  16.06 

4580 

-39.12 

2118 

-15.67 

5280 

-  18.46 

4763 

-40.04 

2122 

-15.49 

5463 

-  19.72 

4025 

-40.48 

2126 

-  17.01 

5500 

-  22.65 

5012 

-41.58 

2130 

19.75 

0 

-20.71 

5012 

-42.47 

2134 

-  20.61 

563 

22.11 

5200 

-42 . 59 

REPRODUCED  AT  GOVERNMENT  EXPENSE 


EPS  :  LOWER  HYDROPHONE 


ANALYIS  METHOD:  QUICK  ANALYSIS 

DATE  OP  ANALYSIS:  U/30/86 


RECORDED  HAGNITODE  LEVELS 


I 


(dB  re  IV) 


OELAYf  MS  )  DIRECT  KACN  DIRECT  PRQ(  UZ  )  REEL  KAON  REEL  PRQ(  UZ  )  NOISE  KAGN 


i 


8 


1950 

19.75 

0 

19.75 

0 

-53.79 

1954 

-12.32 

725 

19.75 

0 

-53.47 

1950 

-11.40 

000 

19.75 

0 

-52.60 

1962 

-12.43 

925 

19.75 

0 

-51.99 

1966 

-12.62 

1063 

19.75 

0 

-50.91 

1970 

-10.87 

1150 

19.75 

0 

-50.04 

1974 

-9.40 

1275 

-19.97 

537 

-48 .96 

1970 

-7.57 

1325 

-20.77 

675 

-40.13 

1902 

-9.42 

1308 

-19.33 

800 

-48.05 

1986 

-13.17 

1500 

-17.37 

0SO 

-49.74 

1990 

-14. 19 

1688 

-17.16 

1000 

-49.13 

1994 

-12.49 

1730 

-16.64 

1075 

-48.17 

1998 

-13.37 

1813 

-16.99 

1150 

-48.32 

2002 

-15.20 

192S 

-19.56 

1275 

-49.24 

2006 

-15.11 

2050 

-20.34 

1363 

-49.33 

2010 

-16.04 

2138 

-23.41 

1475 

-49.32 

2014 

-17.56 

2230 

-21.40 

1600 

-49.20 

2018 

-  18.60 

2375 

-19.95 

1725 

-49.21 

2022 

-19.70 

2425 

-16.20 

1775 

-49.95 

2026 

-  19.58 

2600 

-17.47 

1837 

-49.09 

2030 

-  18.23 

2700 

-21.19 

1950 

-49.96 

2034 

-16.74 

2775 

-26.30 

2025 

-49.72 

2038 

-  15.37 

2875 

-25.22 

2225 

-49.12 

2042 

-  15.73 

2963 

-27.51 

2325 

-49.24 

2046 

-15.33 

3075 

-30.01 

2425 

-49.31 

2050 

-  16.37 

3130 

-26.34 

2487 

-48.66 

2054 

-  17.29 

3  300 

-25.26 

2588 

-48.14 

2058 

- 16 . 97 

3400 

-24.01 

2725 

-47.54 

2062 

-17.63 

3463 

-21.72 

2850 

-47.55 

2066 

-18.95 

3600 

-21.47 

2950 

-47.34 

2070 

-19.51 

3700 

-18.51 

3075 

-46 . 64 

2074 

-19.31 

3  800 

-16.43 

3088 

-46.81 

2078 

-  19.03 

3900 

-20.92 

3180 

-47.17 

2082 

-20.04 

3975 

-18.47 

3325 

-47 . 14 

2086 

-20.06 

4125 

-21.14 

3362 

-47.27 

2090 

-19.76 

4225 

-26.53 

3500 

-47.24 

2094 

-19.12 

4325 

-24.27 

3663 

-46.85 

2098 

-  19.03 

4438 

-25.00 

3763 

-46 .63 

2102 

-19.01 

• 

4525 

-21.41 

3838 

-46.00 

2106 

-18.64 

4625 

-20.66 

3930 

-45 . 49 

2110 

-  19.25 

4713 

-23.11 

4088 

-45 . 19 

2114 

-  19.41 

4813 

-21.31 

4150 

-45.55 

2118 

•  19.40 

4900 

-22.86 

4188 

-46. 10 

2122 

-  19.56 

5025 

-29.27 

4268 

-46 . 43 

2126 

-  18.20 

5162 

-  32.61 

4500 

-46.21 

2130 

-17.17 

5262 

-28.85 

4580 

-46.76 

2134 

-  17.76 

5330 

-28.89 

4662 

-47.02 

2138 

-19.66 

5400 

-29.68 

4750 

-48.27 

2142 

-21.11 

5500 

-27.05 

4913 

-40.34 

2146 

19.75 

0 

-22.19 

5000 

-40.96 

2150 

19.75 

0 

-19.67 

5075 

-48.09 

2154 

19.75 

0 

-19.65 

5100 

-48.86 

2158 

19.75 

0 

-20.04 

5238 

-49.20 

2162 

19.75 

0 

-23.01 

5363 

-49.70 

reproduced  at  government  extense 


NPS 


UPPER  HYDROPHONE 


ANALYtS  METHOD:  NORMAL  ANALYSIS 

DATE  OP  ANALYSIS:  12/1/86 


DEtAY(  MS  ) 

RECORDED 

DIRECT  KfcGN 

MAGNITUDE  UVELS 

DIRECT  PRQ<  FTZ  ) 

( dB  rm 

REEL  KAGN 

IV) 

REEL  rRQ(  HZ ) 

NOISE  MAGN 

1900 

19.75 

0 

19.75 

0 

-47  29 

1916 

19.75 

0 

19.75 

0 

-49.07 

1932 

-9.  14 

775 

19.75 

0 

-46 . 88 

1948 

-10.98 

638 

19.75 

0 

-47 . 79 

1964 

-13.94 

787 

19.75 

0 

-47 . 37 

1900 

-11.66 

800 

19.75 

0 

-46.11 

1996 

-9.63 

963 

19.75 

0 

-46-98 

2012 

19.75 

0 

19.75 

0 

-46.98 

2028 

19.75 

0 

19.75 

0 

-46.98 

2044 

19.75 

0 

19.75 

0 

-50.59 

2060 

-7.61 

1386 

-13.54 

725 

-41  32 

2076 

-10.37 

1738 

-15.72 

1025 

-40.08 

2092 

19.75 

0 

19.75 

0 

-40  08 

2108 

-10.00 

1716 

-15.14 

1063 

-40.72 

2124 

-12.50 

2050 

-12.39 

1325 

-30  68 

2 140 

19.75 

0 

19  75 

0 

-38  68 

21S6 

-12.03 

2088 

-9.50 

1375 

-39  19 

2122 

-13.90 

2350 

-13.87 

1713 

-41 . 26 

2108 

19.75 

0 

19.75 

0 

-41.26 

2204 

-13.13 

2325 

-14.15 

1688 

-41  00 

2220 

-16.85 

2613 

-13.45 

1913 

-39.81 

2236 

-14.03 

2438 

-14.45 

1775 

-41  69 

2252 

15.76 

2787 

-11.38 

1975 

-39.08 

2268 

-  13.21 

2975 

-19.75 

2225 

-42  37 

2284 

19.75 

0 

19.75 

0 

-42  37 

2300 

19.75 

0 

19.75 

0 

-42 . 37 

2316 

19.75 

0 

19.75 

0 

-42 . 37 

2332 

-15.87 

3275 

-22.67 

2588 

-40.30 

2  348 

-17.11 

3500 

-20.52 

2825 

-42.85 

2364 

•  15.06 

3362 

-20.24 

2662 

-41  32 

2380 

-17.99 

3650 

-17.01 

2988 

-39  05 

2396 

19.75 

0 

19.75 

0 

\ 

u 

« 

O 

V* 

2412 
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